Analytic Determination of the Critical Coupling for Oscillators in a Ring 
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We study a model of coupled oscillators with bidirectional first nearest neighbours coupling with 
periodic boundary conditions. We show that a stable phase-locked solution is decided by the oscil- 
lators at the borders between the major clusters, which merge to form a larger one of all oscillators 
at the stage of complete synchronization. We are able to locate these four oscillators as well as the 
size of major clusters in the vicinity of the stage of full synchronization which we show to depend 
only on the set of initial frequencies. Using the method presented here, we are able to obtain an 
analytic form of the critical coupling, at which the complete synchronization state occurs. 
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Weakly coupled oscillators in the chaotic state 
have been known to represent many physical sys- 
tems, as well as chemical, biological, neurologi- 
cal and so on. These systems synchronize in fre- 
quency under the influence of coupling. Knowing 
beforehand the value of the coupling constant and 
the dynamical behavior of the individual oscilla- 
tors for complete synchronization to occur is an 
important source of information for real applica- 
tions. This paper is a continuation of previous 
theoretical results for these systems. Here, we 
derive relationships that allow us to determine 
the oscillators which first lock in phase and drag 
the whole system into the synchronized state as 
well as the size of the two existing clusters before 
the transition. 



I. INTRODUCTION 

In recent years we have seen oscillators coupled 
through nearest neighbors interactions to be used to un- 
derstand the behavior of systems in physics, chemistry, 
biology, neurology as well as other disciplines, to model 
several phenomena such as: Josephson junction arrays, 
multimode lasers, vortex dynamics in fluids, biological 
information processes, neurodynamics 1 -^' 3 -! 4 -. These sys- 
tems have been observed to synchronize themselves to a 
common frequency, when the coupling strength between 
these oscillators is increased 4 *^. In spite of the diversity 
of the dynamics, the synchronization features of many of 
the above mentioned systems might be described using 
a simple model of weakly coupled phase oscillators such 
as the Kuramoto model 4 ^^, as well as its variations to 
adapt it for finite range interactions which are more re- 
alistic to mimic many physical systems. Difficulties arise 
since finite range coupled systems are difficult to ana- 
lyze and to solve analytically. In spite of that, in order 
to figure out the collective phenomena when finite range 



interactions are considered, it is of fundamental impor- 
tance to study and to understand the nearest neighbour 
interactions, which is the simplest form of the local inter- 
actions. In this context, a simplified version of the Ku- 
ramoto model with nearest neighbour coupling in a ring 
topology, which we shall refer to as locally coupled Ku- 
ramoto model (LCKM), is a good candidate to describe 
the dynamics of coupled systems with local interactions. 
Several reports exist where the LCKM has been used to 
represent the dynamics of a variety of systems, such as 
Josephson junctions, coupled lasers, neurons, chains with 
disorders, multi-cellular systems in biology and in com- 
munication systems^^. It has also been shown that 
the equations of the resistively shunted junction which 
describe a ladder array of overdamped, critical-current 
disordered Josephson junctions that are current biased 
along the rungs of the ladder can be expressed by a 
LCKMi£. For nearest neighbours coupled Rossler os- 
cillators the phase synchronization can be described by 
the LCKMii, as well as locally coupled laser a 12 i 13 , where 
local interactions are dominant. LCKM can also be used 
to model the occurrence of travelling waves in neurons 4 ^. 
In communication systems, unidirectionally coupled Ku- 
ramoto model can be used to describe an antenna ar- 
rayi 4 -. Such unidirectionally coupled Kuramoto models 
can be considered special case of the LCKM and 
it often mimics the same behaviour. Therefore, LCKM 
can provide a way to understand phase synchronization 
in coupled systems in general. 

While in the Kuramoto model for long range interac- 
tions one has to rely on average quantities, in a mean 
field approximation or by means of an order parameter, 
etc., in the local model it is necessary to study the be- 
haviour of individual oscillators in order to understand 
the collective dynamics. Therefore, due to the difficulty 
in applying standard techniques of statistical mechanics, 
one should look for a simple approach to understand the 
coupled system with local interactions by means of nu- 
merical study of a temporal behaviour of the individual 
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oscillators. Such analysis is necessary in order to ob- 
tain a close picture of the effect of the local interactions 
on synchronization. In this case, numerical investiga- 
tions can assist to figure out the mechanism of interac- 
tions at the stage of complete synchronization which in 
turn helps to get an analytic solution. Earlier studies 
on the LCKM show several interesting features including 
tree structures with synchronized clusters, phase slips, 
bursting behaviour and saddle node bifurcation and so 
o n 15 ' 16 . It has also been shown that neighbouring ele- 
ments share dominating frequencies in their time spectra, 
and that this feature plays an important role in the dy- 
namics of formation of clusters in the local modelii; that 
the order parameter, which measures the evolution of 
the phases of the nearest neighbour oscillators, becomes 
maximum at the partial synchronization points inside the 
tree of synchronization^ and a scheme has been devel- 
oped based on the method of Lagrange multipliers to 
estimate the critical coupling strength for complete syn- 
chronization in the local Kuramoto model with different 
boundary conditions^. 

Very recently, we identified two oscillators which 
are responsible for dragging the system into full 
synchronization^, and the difference in phase for this 
pair is ±7r/2. In this work we develop a method to obtain 
an analytic solution for the value of the critical coupling 
at which full synchronization occurs, once a set of initial 
conditions for the frequencies of the N oscillators is as- 
signed. This method will allow us not only to calculate 
the analytic form of the critical coupling but also to de- 
termine the number of oscillators at the major clusters 
in the vicinity of the critical coupling as well as to de- 
termine which is the pair of oscillators that has a phase 
difference ±7r/2 at the stage of full synchronization. 

This paper is organized as follows. In Sec.[TI]we investi- 
gate the LCKM where periodic boundary conditions are 
used. We derive an analytic form for the critical coupling 
at the stage of complete synchronization as well and de- 
termine the number of oscillators at each cluster in the 
vicinity of the critical coupling. Finally, in Sec. IIIII we 
give a conclusion which is based on a summary of the 
results. 



II. OSCILLATORS IN A RING 

The local model of nearest neighbour interactions, or 
LCKM can be considered as a diffusive version of the 
Kuramoto model, and it is expressed a a 16 ' 17 i 18 ' 19 ' 20 : 

9i = u>i + A"[sin(0j) - sin(0<_i)], (1) 

with periodic boundary conditions 8i+N = Qi and for 
2 = 1,2, N. The set of the initial values of frequencies 
{u)i} are the natural frequencies which are taken from a 
Gaussian distribution and K is the coupling strength. 
The phase difference is defined as fa = 6 i+ i — 9i for 
i = 1, 2, 3, N. These nonidentical oscillators |T]) clus- 
ter in time averaged frequency, until they completely syn- 



chronize to a common value given by the average fre- 
quency luq = jj^Z^i^i, at a critical coupling K c . At 
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FIG. 1: Synchronization tree for a system of 30 oscillators, 
with detailed composition of each cluster before full synchro- 
nization. 

K > K c the phases and the frequencies are time inde- 
pendent and all the oscillators remain synchronized. In 
Fig.l we show the synchronization tree for a periodic sys- 
tem with N = 30 oscillators, where the elements which 
compose each one of the major clusters are indicated in 
each branch. These clusters merge into one at K c where 
all oscillators have the same frequency. The major clus- 
ters just before K c contain Ni and N2 oscillators, where 
N = Ni + N 2 . It is not necessary for these clusters to 
have the same numbers of oscillators. At the vicinity 
of K c , major clusters of successive oscillators have sets 
of nearest neighbours at the borders. An interesting fact 
emerges: the phase-locked solution is always valid for one 
and only one phase difference, and this phase difference 
is between two oscillators at the border of the clusters^. 
Thus, for these two neighbouring oscillators, the equation 
for the phase difference is: 

4> n = A„ - 2Ksin(0 n ) + K sin(^„_i) + K sin(0„+i), 

(2) 

where A„ = uo n +i — w n . Equation ^ at K c has <fi n = 0, 
and hence 9 n = 9 n+ i = ujq. It has been found that 
the phase-locked solution is satisfied when <j> n = ir/2 
for the case of w n +i > u) n and 4> n = — tt/2 for the re- 
verse. In addition the phase-locked solution exists al^ 
X n = \j^ + [sin(0 n _i) +sin(0 n _i)] | = 2. It is already 
well known that in the vicinity of K c , phases as well as 
frequencies present the phenomenon of phase slip, that 
is, they remain constant for a given period of time T and 
then they jump, followed by a another period of con- 
stant value, a jump and so on. During this period T 
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equation (J2J) can be solved analytically to give quanti- 
ties (f> n and (j) n . Due to the diffusive character of the 
LCKM, all other quantities ifi's and 4>'s of other oscilla- 
tors relate to </>„ and </>„ (see details in reference^), and 
will in turn present the same phenomena, which has also 
been seen numericall y 11 ^ 6 . Thus, if one can develop a 
method to allocate the two oscillators which will have 
the phase difference <j> n = ±7r/2, then it will be possible 
to determine the value X n and to obtain the critical cou- 
pling at which a complete synchronization occurs (K c ). 
A difficulty arises due to the determination of the values 
of sin(</>„ + i) and sin(</>„_i) in addition to the topology 
of a ring which is an endless system. Therefore, there 
is no direct method to specify the two oscillators n and 
n + 1, which have phase difference that would satisfy the 
phase-lock condition. 
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FIG. 2: Selected Values of sin <j>i at K = K c for the system of 
30 oscillators of Fig.l. 

We can take advantage that there are four oscilla- 
tors, now labeled I, I + 1, m and m + 1 , at the borders 
of the major clusters in the vicinity of K c from which 
only one pair will have a phase difference correspond- 
ing to the phase-locked solution; i.e, |sin(</>„)| = I. As 
shown in Fig. 2, the values of sin(</>/) and sin((/) m ) are 
always the maximum and minimum of the sin(<^)'s for 
all phase differences. From these two phase differences 
4>i and <p m , one of them has a value ±7r/2, while the 
other would be close T 7r /2, getting closer as N increases. 
This fact has been verified numerically for several real- 
izations of N and {w^j's. If we start adding equations 
of the systems in a ring JT]) (adding elements) in a simi- 
lar way as in reference^, we generate a sequence Zi , with 
Zi = iu}Q — Y^j=x f° r * = 2, 3, N — 1. After a detailed 
study of existing correlations, we arrive to a criteria to 
determine the four oscillators at the borders of the major 
clusters in the vicinity of K c . We find that the maximum 
value of Zi refers always to the oscillators at one border 
while the minimum of Zi points to the oscillators at the 
other border of the given cluster. One of them is always 
positive and the other is negative. The sign of Zi de- 
pends on the values of Aj = lui+i — ujf. for max(Zi), 



Aj > 0, while min(Zi) corresponds to Aj < 0. A thor- 
ough study also shows that the value of sin(0j) > cor- 
responds to max(Zi) while the sin(^) < corresponds 
to min(Zi). Thus we can calculate all the Z[s and assign 
the maximum and minimum values which refer to the in- 
tegers / and m. However, we have not yet resolved which 
one refers to the two oscillators with the phase difference 
equals ±7r/2. 

Since now we know how to point to the four oscillators 
I, I + I, m and m + 1 at the borders of major clusters 
at K Cl it is possible to obtain an analytic form for K c , 
within reasonable accuracy. These two quantities Zi and 
Z m are related to both values of 4>i and <p m for a system 
of coupled oscillators in a ring ([I]) such that: 



K c [sin(0j) - sin(0jv)] = Z% 
K c [sin(0 m ) - sm(<f) N )] = Z m . 



(3) 
(4) 



The value of sin(0jv) presents a difficulty in determining 
which one of the phase differences (pi or <j) m corresponds 
to \tt/2\, since the phase- lock condition at the critical 
coupling will bei£ either K c = Zi + K c sm((f>N) or K c = 
Z m +K c sm(4>N). We can then rely on numerical findings 
and characteristics of 4>i and <j) m to use equations (3) 
and (4) to obtain an approximate analytic expression for 
K c . A detailed numerical investigation of sin(0jv) shows 
that it is always small, in comparison to both values Zi 
and Z m , and that K c depends mainly only on these two 
quantities Z; and Z m . Taking this fact into consideration 
and since both sin(^) and sin(0 m ) are always opposite 
in sign as well as the two quantities Zi and Z m , we use 
equations (3) and (4) to obtain 



K, 



\Zt\ 
2 



(5) 



where e depends on the difference between both | sin(0/)| 
and |sin(0 m )|. Helping ourselves by numerical studies 
we find that e depends on the quantities A/ and A m and 
it can can be written as e w |A||+|A,„|-|A,+A m ^ ^phus we 

16 

obtain an approximate expression for K c , which we call 
K£, and is given by: 



\Zi\ + \Z„ 



I A; + A I; 
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(6) 



Fig. 3 summarizes numerical investigations of the deter- 
minations of K c . We plot log-Kc versus log A from nu- 
merical simulations of equation ([T]) (triangles) and from 
the analytic results given by ([6]), first by considering e 
going to zero (squares) and then taking its complete de- 
pendence on both A; and A m (circles). The validity of 
equation (6) is clearly shown for values of iV ranging from 
30 to 1000. The dependence of K c on both \Z t \ and \Z m \ 
as in equation © and as N increases becomes clear. It 
can also be inferred that the term which depends on e 
becomes negligible. This is due to the fact that as N 
increases the oscillators of indexes I and I + 1 are be- 
coming closer in frequencies to each other as well as the 
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FIG. 3: (Color online) log-log plot of K c versus N from nu- 
merical simulation of system (1) (triangles) and analytic cal- 
culation using equation (6) (circles) and equation (5) with 
e = (squares). 



two oscillators of indexes m and m + 1 . We also observe 
that K c grows as ~ O^J~N , in the same limit as found 
by Strogatz and Mirolloi^ (for details, sec explanation in 
this reference). 

Summarizing, if one knows the set of initial frequencies 
{tOi}, it is possible to point at the four oscillators at the 
borders of the major clusters just below K c and then the 
calculation of K c is performed using equation (6) (thus 
obtaining K®), without the need of computer simulation 
of system ([T]), just using the values of Z\ and Z m . If we 
are interested in determining which phase difference will 
have a phase-lock condition ±7r/2, we use the fact that 
sin(</>/) and sin(<^ m ) have opposite signs as well as they 
are maximum and minimum among all values of sine of 
the phase differences. The sign of the quantity sin(^)jv) 
has the same sign of the quantity X\ = — {Zi + Z m )/2, 
which is taken from the sum of equations (3) and (4) 
(eliminating for a moment the small difference between 
sin(0/) and sin((/> m )). Depending on the signs of Z\ and 
Z rn , we know the signs of sin(<^) and sin(</> rra ), and hence 
the sign of sin(<^/v)- Therefore, we count two quantities 
x 2 = ±K c c l — Z\ and x 3 = ±K® — Z m , positive sign for 
Z > and negative sign for the reverse. Three cases 
will exist: first from the quantities x 2 and £3, one is 
positive and the other is negative. Thus depending on 
the sign of x\ we choose either x 2 or X3 to be K® sin(aW). 
Second x 2 and X3 have the same signs, then we check the 
minimum between \x\ — x 2 \ and \x\ — Xa\ and depending 
on which one is the minimum, we take either x 2 or X3 
to be if"sin(0jv)- Third \x\ — x 2 \ = \x\ — 2:3 1, then 
we take the minimum outcome of x 2 and £3. Now we 
know the value of if°sin(^>jv) and its sign. Therefore, 
we know which equation from (3) and (4) will be used to 
give K" . Thus we specify which phase difference of index 



/ or m would have ±7r/2. We tested this method on the 
simulations we have done and it matches the outcome of 
the numerical simulations. 

The number of oscillators in each cluster at the vicinity 
of K c can be determined, once we assigned the indexes 
I, I + 1, to and m + 1, which, we remind the reader, 
are obtained from Z[ and Z m , maximum and minimum 
values of the sequence Zi. The size of one cluster of 
Ni oscillators is determined by counting the difference 
Ni = (to + 1) — I and the size of the other cluster is de- 
termined as N 2 — N — Ni . Similarly to the calculation of 
Zi and Z m7 we can determine other two quantities which 

are Yi = Nxuj - w i and Y 2 = N 2 u - X)i=m+i LJi 

taking into consideration the periodic boundary condi- 
tions. It is found that |Yi| = \Y 2 \. These quantities are 
related to Zi and Z m by Y\ = Z m —Z\ = —Y 2 . It is easy to 
show that K c = + e = + e. The two quantities Y\ 
and Y 2 provide a criteria to understand synchronization- 
desynchronization at K c . If one arrives from above K c 
where all oscillators are synchronized and have the same 
value of frequency, at K c the oscillators split into two 
groups of N\ and N 2 , at K c , depending on these two 
quantities Y\ and Y 2 , where \Yi\ = \Y 2 \. It is not neces- 
sary for N± to be equal to N 2 . Both quantities Y\ and Y 2 
have opposite signs since they refer to two groups of os- 
cillators (two clusters) one of them rotates with average 
frequency over loq and the other has an average frequency 
lower than ujq. 

Comparing our findings of K c with the work of 
Daniels et al.— , our method has the advantage of finding 
the value of K c without performing numerical simula- 
tions once we know the set of initial frequencies {w;}. 
In addition we get the condition of synchronization- 
desynchronization at K c and obtain the number of os- 
cillators in each branch in the vicinity of K c . 



III. CONCLUSION 

We have analyzed the conditions of the phase differ- 
ences for the onset of complete synchronization at the 
critical coupling strength in a Kuramoto-like model with 
nearest neighbour coupling with periodic boundary con- 
ditions. Such analysis allows us to determine the four 
oscillators located at the borders of the major clusters 
(formed by successive oscillator) which will meet at the 
critical coupling to form one cluster of all synchronized 
oscillators. With the help of these findings we derive 
the analytic expression for the critical coupling when all 
oscillators will have the same frequency and, phase dif- 
ferences and instantaneous velocities become time inde- 
pendent. In addition, we are able to determine which is 
the phase difference, that will have a phase-lock solution 
±7r/2. From the derivation we also extract the size of 
the clusters before complete synchronization. The ana- 
lytic form of K c depends only on the initial frequencies, 
through the quantities Zi and Z m , where the indexes I 
and to correspond to the borders of the clusters. The 
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quantities Zi and Z m correspond to the maximum and 
minimum values of the sequence Z^. These quantities 
in fact are related to the statistics of the distribution 
of the set of initial frequencies {coi}, when this sample is 
obtained from a Gaussian distribution, as shown by Stro- 
gatz and Mirollo^. A detailed study within this context 
could shine light on the behavior of K c as N — > oo, not 
just for the case of a Gaussian distribution, but for others. 
This investigation plus extension of the method to study 
cluster formation inside the tree will be topics of further 
analysis. The advantages of the study presented here 
is that we can determine the value of the coupling con- 
stant that will synchronize the system of coupled oscil- 



lators without carrying out numerical simulation as well 
as to determine the sizes of the clusters just before this 
happens. Generalization of these results to different cou- 
plings and boundaries are under investigation and will be 
presented elsewhere. 
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